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SUMMARY 


A program to investigate the aerodynamic heat transfer of a non- 
isothermal hemisphere-cylinder has been conducted in the Langley 11-inch 
hypersonic tunnel at a Mach number of 6.8 and a Reynolds number from 
approximately 1.09 X 10^ to 1.03 X 10^ based on diameter and free-stream 
conditions. The experimental heat-transfer coefficients were slightly less 
over the whole body than those predicted by the theory of Stine and Wanlass 
(NACA Technical Note 33^0 for an isothermal surface. For stations within 
45° of the stagnation point the heat-transfer coefficients could be cor- 
related by a single relation between local Stanton number and local 
Reynolds number. 

Pitot profiles taken at a Mach number of 6.8 on a hemisphere cylin- 
der have verified that the local Mach number or velocity outside the 
boundary layer required in the theories may be computed from the surface 
pressures by using isentropic flow relations and conditions immediately 
behind a normal shock. The experimental pressure distribution at a Mach 
number of 6.8 is closely predicted by the modified Newtonian theory. The 
velocity gradients calculated by using modified Newtonian theory at the 
stagnation point vary with Mach number and are in good agreement with 
those obtained from measured pressures for Mach numbers from 1.2 to 6.8. 

At the stagnation point, the theory of Sibulkin, using the diameter 
and conditions behind the normal shock, was in good agreement with the 
experiment when the velocity gradient at the stagnation point appropriate 
to the free-stream Mach number was used. 


INTRODUCTION 


One of the earliest experimental investigations in connection with 
heat transfer on a hemisphere-cylinder at supersonic speeds was carried 
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out by Eber (ref. l) who obtained the "brake" temperature distribution 
along the surface. These brake or effective recovery temperatures 
decreased with increasing distance from the stagnation point on the 
hemisphere and were influenced by the free-stream Mach number . Korobkin 
(ref. 2) and Stine and Wanlass (ref. 3) show the same type of distribu- 
tions and, in addition, demonstrate that for the hemisphere-cylinder 
with a laminar boundary layer the recovery factor t) may be expressed 

by the square root of the Prandtl number based on the local conditions 
just outside the boundary layer. 

Stalder and Nielsen (ref. 4) compared the average heat transfer 
(Nusselt number) with Reynolds number for different Mach numbers and 
found that, in their experiments, data at subsonic and low supersonic 
Mach numbers could be correlated when the supersonic data were based on 
conditions behind the normal shock. Experimental investigations of the 
local heat transfer on a hemisphere-cylinder have been made by Korobkin 
(refs. 2 and 5) and Stine and Wanlass (ref. 3) • Stine and Wanlass 
noticed that, as the free-stream Mach number was increased, the pressure 
distributions were tending toward a common curve and suggested that the 
heat-transfer results they obtained at a Mach number of 1.97 could be 
representative of those obtained on similar isothermal surfaces with 
laminar boundary layer for Mach numbers greater than 1.97 and of tempera- 
tures less than that of dissociation. Although this tendency is correct, 
a Mach number of 1.97 is probably too low. 

Some theoretical methods for obtaining local heat transfer on a 
hemisphere-cylinder at supersonic speeds have as their bases the solu- 
tions of incompressible flow on isothermal surfaces. An exact incom- 
pressible solution utilizing a method applied by Squire (presented in 
ref. 6) to a cylinder was given by Sibulkin in reference 7 for the heat 
transfer at the stagnation point. Sibulkin suggested that this solution 
may be applicable at supersonic speeds when conditions behind the normal 
shock are used. An approximate solution was derived for incompressible 
flow over a hemisphere with an isothermal surface by Sibulkin (presented 
in ref. 2) by assuming that the velocity and temperature profiles were 
similar and of parabolic shape . Korobkin (ref. 2) compared these theo- 
retical results with his data at a Mach number of 2.80 and found that, 
when the incompressible-flow solution was adjusted so as to pass through 
Sibulkin’ s stagnation-point value, it roughly approximated the data. 
However, this result does not agree with the theory by Stine and Wanlass 
(ref. 3) who used a different approach. By using the Mangier transforma- 
tion, the Stewartson transformation, and thermal solutions to the Falkner- 
Skan wedge-flow problem, they evaluated the heat-transfer rate in axisym- 
metric flow for an isothermal surface in terms of the known heat-transfer 
rate in an approximately equivalent two-dimensional flow. 

The purpose of this investigation was to extend the range of the 
previous investigations for heat transfer on hemisphere-cylinders to a 
Mach number of 6.8. This investigation utilized data from transient 
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noni sot hernial temperature distributions since, at the high temperature 


levels 

involved, an isothermal surface is difficult to obtain. 

SYMBOLS 

a 

speed of sound 

b 

semihemi spherical arc length 

c m 

specific heat of model wall material 

c p 

specific heat of air at constant pressure 

c v 

specific heat of air at constant volume 

D 

diameter 

h 

local aerodynamic heat-transfer coefficient, = — ^ - - 

l r - I s 

k 

thermal conductivity for air 

V 

thermal conductivity for model 

L 

a characteristic length 

M 

Mach number, u/a 

Nu 

Nusselt number, hL/k 

Pr 

Prandtl number, c p p/k 

P 

pressure 

q. 

heat flow to model per unit area, -k — ) 

Ws 

Re 

T ' * 

Reynolds number, 

P 

r 

radius 

St 

Stanton number based on x and local free-stream conditions 
just outside boundary layer, St = Nu/RePr = h/puc p 


h 
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T 

absolute temperature 

u 

veloc ity 

X 

distance along body generator from stagnation point 

y 

space coordinate normal to body 

Z 

distance along cylinder from juncture 

P 

velocity gradient along surface, ^ 

e 

emissivity 


T r - T-l 

recovery factor, 

*t “ T 1 

9 

angle from stagnation point on sphere 

M- 

absolute viscosity 

V 

kinematic viscosity 

1 

angle about axis of model 

P 

mass density of air 

Pm 

mass density of model shell material 

a 

Stephen-Boltzmann constant 

T 

time 

Subscripts: 


a 

condition of model's surroundings 

D 

parameter based on static conditions immediately behind normal 
shock and with model diameter as a characteristic length 

h 

pitot tube behind bow shock 

i 

not corrected for heat conduction along model skin 

in 

local condition on inside surface of model shell 

L 

parameter using L as characteristic length 

r 

recovery, surface temperature 
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s local conditions on outside surface of model 

t stagnation conditions 

x parameter based on static conditions just outside 

boundary layer and with x as characteristic length 

a static conditions immediately behind normal shock 

1 local conditions just outside boundary layer 

00 undisturbed free-stream conditions 


APPARATUS 

Tunnel 


The present investigation was conducted in the Langley 11-inch 
hypersonic tunnel described in references 8 and 9* Air is stored at 
50 atmospheres pressure and is released through an adjustable pressure 
regulating valve and a new direct air heater with tube resistance elements 
of nickel-chromium alloy (replacing the storage heater described in 
refs. 8 and 9) to the settling chamber and nozzle. A two-dimensional 
nozzle constructed of Invar and designed for a Mach number of 7 was used 
for this investigation. Invar was used in the construction of this 
nozzle in order to alleviate the deflection of the first minimum, which 
occurred in the steel nozzle because of differential heating of the 
nozzle blocks. Preliminary calibrations have indicated a Mach number 
of 6.86 ± 0.04- at 30 atmospheres stagnation pressure in the central core 

of uniform flow which measures about 6~- inches in the vertical direction 

by about 6 inches in the horizontal direction. 


Tunnel Conditions 

During the tests the tunnel was operated through a stagnation- 
temperature range from 1,04-0° R to l,l60° R and through a stagnation- 
pressure range from 10 to 31 atmospheres. The free-stream Reynolds 
numbers (based on body diameter) were in the range from approximately 

3.4- X 105 to 1.03 x 10^. The calibrated tunnel Mach number at the model 
nose was used to correlate the data. In these tests, the model nose was 
approximately 6 inches ahead of the center of the test section. At this 
point, the tunnel calibration gives a Mach number of 6.82 for a stagna- 
tion pressure of 31 atmospheres, which agrees with values obtained from 
pressures at the nose of the pressure model. As the pressure decreases. 
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the Mach number also decreases slightly, so that, at 10 atmospheres stag- 
nation pressure, the Mach number 6 inches ahead of the test-section 
center line is 6 . 7 ^-. 

Instrumentation 

The temperature measurements included stagnation temperature and 
model skin temperature. The stagnation temperature was measured by 
shielded thermocouples distributed at various stations in the settling 
chamber . 

The potentiometers used for recording the model temperature showing 
greatest transient rate of change were cyclic printers with a 1-second 
interval, and with a 1-second full-scale response. In order to reduce 
uncertainties associated with the time response of the recording instru- 
ments, the thermocouples on the sphere were connected to be read in order 
of decreasing temperature about the sphere from station 1 to station 11. 
The temperature at station 1 was recorded twice to allow time for bal- 
ancing out the large temperature change between stations 1 and 11. Four 
instruments were connected to read the temperatures of the first 11 ther- 
mocouples in this manner and were synchronized so that each thermocouple 
was read once each 3 seconds by one of the potentiometers. 

Model surface pressures were recorded on film by the evacuated 
capsule instruments described in reference 8. The motion of a diaphragm 
rotates a small mirror to displace the trace of a light beam falling on 
a moving film. Pressure cells were chosen to give as near full-scale 
deflection as possible for the measuring station. The stagnation pres- 
sure was measured with high-accuracy Bourdon tube gages. 


Schlieren System 

Two schlieren systems were used in this investigation. One system 
had a single-pass, vertical Z-light path with a horizontal knife edge, 
and the other, a double-pass horizontal light path with the light source 
and cut-off knife edge at the center of curvature of the concave spherical 
mirror. The double-pass system was used only when greater sensitivity 
was needed for boundary -layer study. Film exposures were of approximately 
3 microseconds duration. 


Models 


The four models used in this investigation are shown in figures 1 
and 2. The pressure-distribution and heat-transfer models were constructed 
of SAE 1020 carbon steel and the temperature-recovery model was constructed 
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by electroplating nickel on a mandrel. The dimensions of the heat- 
transfer models were carefully measured before the parts were assembled 
with silver solder. 

The large temperature-distribution model was 3.025 inches in diam- 
eter, 10.5 inches in length, and varied from O.O98 inch to 0.101 inch in 
wall thickness. It contained 25 thermocouples, silver soldered into holes 
in the model skin so that their effective measuring junctions were located 
on the inner surface. These thermocouples were of Wo. 3 6 gage chromel 
alumel wire and were arranged in a spiral of one revolution on the 
hemisphere and of two revolutions on the cylinder. Ten of these ther- 
mocouples were located on the hemisphere, fourteen were on the cylinder, 
and one was on the sphere -cylinder juncture as given in figure 1. 

The pressure-distribution model was 3.000 inches in diameter, 

10.5 inches in length, and 0.10 inch in wall thickness and contained 
16 pressure orifices located as shown in figure 1. Two of the orifices 
were located on the base of the model. The orifices were 0.0^0 inch in 
diameter and were arranged in a spiral of one revolution about the sphere 
in order to reduce any orifice-interference effect. The orifices on the 
cylinder were also arranged in a spiral of one revolution. 

The small temperature-distribution model was I.I70 inches in diam- 
eter, 5.32 inches in length, and varied from 0.030 inch to 0.035 inch in 
wall thickness. Temperature measurements on this model were taken from 
eight chromel alumel thermocouples of No. 30 gage wire, five of these 
being located at stations from the stagnation point to the sphere -cylinder 
juncture. 

The temperature -recovery model constructed of electroplated nickel 
was 2.000 inches in diameter, 5. 00 inches in length, and 0.005 inch thick. 
Temperature measurements were taken from it chromel alumel thermocouples 
of No. 30 gage wire. 


ANALYSIS 

Heat Transfer to Model 


The local rate of heat flow into the model was calculated from the 
local temperatures as a function of time and position. The general heat- 
conduction equation (ref. 10) in a homogeneous material is 

ST = V 2 t 

St P m c m 


( 1 ) 
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In spherical coordinates this equation becomes: 


dT = 


1 ^ 1 

>2 dT\ 

+ 1 .a 

/sin 0 

1 d 2 T 

St 

p m c m 

r 2 dr' 


r 2 sin 0 

pill O T 

l 

r 2 sin 0 d| 2 


Since the model shell is relatively thin with respect to its other 
dimensions, it is assumed that the first and second derivatives of T 
with respect to 0 and £, as well as the first derivative of T with 
respect to t, are independent of r. Equation (2) may then be integrated 

St 

once with respect to r, using the boundary condition = 0 at the 

dT 

inner surface where r = r. , and then may be solved for “ to obtain 

m dr 

the following: 

l r3 -,,^to£Y-m2m dT _ q St\ _ ~ r in\ dfg ( 3 ) 

dr l 3r 2 j V yr sin 0/d0 \ d0/ \r 2 sin 0/ d| 2 


Equation (3) may then be evaluated at the outer surface where r = r s 

and at zero angle of attack so that T is independent of |j this 
evaluation gives the following: 





| Pm c m dT 

km 8 t 


'in\ d 


[r^sin 0/d0 



However at the surface, 



Therefore , 


ds 


= -p c 
m m 




( 4 ) 


( 5 ) 


( 6 ) 


An examination of the above equation discloses that the first term 
on the right represents q i? the local net heat storage which would be 
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realized in a shell element if the conduction were zero. The term 

.5 - r- 3\ 

2 j is approximately ^r s - r^j , or the shell thickness, so 

that in the limit as t — > 0, the q^ term becomes the same as that used 

for a flat plate. The second term on the right represents the local net 
heat conducted along the model shell. 

When this equation was used to calculate q g at the various times 

dT 

and stations, the values of — for the heat -storage term were taken 

or 

from faired curves of temperature plotted against time. The values of 
J^Ysin 6 ^ j for the conduction term were found graphically by taking 

slopes of faired curves. 


In order to give a comparison of the relative size of the heat- 
storage and conduction terms for the tests using the 3 . 025 -inch- diameter 
heat-transfer model, typical q s and q.j_ values are plotted in figure 3 

for a time early in the run and for a time late in the run. The difference 
between these curves q s - q i is the value of the conduction term. It 

is seen that, at 5 seconds from the start of the run, the conduction term 
is only 20 percent of qg for the station at the stagnation point; but, 
near the end of the rim, the conduction term is about 70 percent of q s 
at the same station. 


The local heat-transfer coefficients were determined from the standard 
heat-flow equation for convection, qg = h ^T s - T r ) . The quantities T r 

and h were assumed constant for a given set of test conditions, and h 
was determined from the slope of the straight-line variation of q g 

plotted against T s as is shown in figure 4. In the present tests, 
there were some variations in the test conditions during a run. The 
stagnation temperature sometimes varied 20°, and the pressure varied 
nearly 5 percent in the lowest pressure run. In order to correct 
partially for the variation in stagnation temperature, the values of 
q s /T t were faired through a calculated value of T t . The recovery 

temperature ratios used in the fairing for the determination of h were 
obtained from the recovery factors for a flat plate given by Van Driest 
(ref. 11) by using the local flow conditions at the outside edge of the 
boundary layer. 


Errors 

The present heat-transfer tests were affected by four small sources 
of error which were of uncertain magnitude and were not directly 
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incorporated in the computations. One source was the heat loss from the 
model by radiation, and another was the heat loss by conduction down the 
thermocouple wire. The third was the assumption that the difference in 
temperature distribution between the outside and the inside surfaces of 
the model was negligible. The fourth source of error is the assumption 
that the heat-transfer coefficient is independent of the temperature 
level. The error associated with radiation was estimated to be about 
3 percent of the value of q at the highest temperature condition. The 
error associated with heat flow down the thermocouple wire was not esti- 
mated directly but was assumed to be small, especially on the large 
heat-transfer model where the thickness of the model was 0.10 inch and 
the diameter of the thermocouple wire was 0.005 inch. The temperature 
difference between the outside and the inside skin surface was less than 
10° F at the stagnation point immediately after the start of a run at 
p_l_ = 31 atmospheres. This temperature difference represented an error 

of less than 0.01 in the determination of T s / T^ . It was noticed, how- 
ever, that each of the first three errors tended to reduce the calculated 
value of the heat -transfer coefficient, so that the final results of this 
data analysis might tend to be slightly low when compared with theory. 

For the assumption that the heat-transfer coefficient is independent of 
the temperature level, the maximum variation in wall temperature occurred 
at the stagnation point and here the wall temperature varied from 180° F 
to 55 0° F. By the theory of Reshotko and Cohen (ref. 12) a variation 
in wall temperature of this magnitude would result in a deviation of 

only tl^j- percent from the mean value of Nu/ N / _ Re at the stagnation point. 


RESULTS AND DISCUSSION 
Pressure Distribution 


The longitudinal pressure distributions for the hemisphere-cylinder 
at zero lift are presented in figure 5 for Reynolds numbers of 0.71 X 105 
and 1.42 X 10^ based on diameter and conditions behind the normal shock. 
These pressure distributions are presented as a ratio of local pressure 
coefficient to the pressure coefficient at the stagnation point. This 
pressure distribution has been compared with the Newtonian theory from 
reference 13, the data of Stine and Wanlass (ref. 3) at M = 1.97 and 
3.80, and the data of Oliver (ref. 14-) at M = 5 *8* It will be noted 
that, at Mach numbers above 3*80, the pressure-coefficient distribution 
curves asymptotically approach a single curve, and the local Mach number 
distribution on the hemisphere becomes independent of the free-stream 
Mach number. This was pointed out in reference 3 and is in accordance 
with the hypersonic similarity parameter (ref. 15 ) which predicts that 
the flow pattern about an arbitrary body will become stationary as Mach 
number increases; for a blunt body, this tendency occurs at a lower Mach 
number than that for a slender body. 
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In this report, the isentropic-flow relations between the pressures 
measured at the stagnation point and the local surface pressures measured 
at subsequent stations were used to calculate local flow conditions just 
outside the boundary layer. The validity of this method was investigated 
by the use of pitot surveys of the boundary layer. In figure 6, the 
pitot profiles are presented for the hemisphere -cylinder juncture, a 
station midway along the cylinder, and a station l/8 inch in front of the 
base of the cylinder. These profiles are characteristic of laminar 
profiles and indicate that the boundary layer remains laminar over the 
whole body. The point where the pitot profiles break away from the 
characteristic boundary -layer profile corresponds closely to the pitot 
ratio obtained by using the static surface pressure from the pressure 
model and the isentropic-flow relations (ref. 1 6) from the stagnation 
point to the survey stations. Thus, the pressures from the body surfaces 
are sufficient through the use of isentropic-flow relations to obtain 
local flow conditions outside the boundary layer for station at least up 

to = k.S. In figures 7 attd 8 , the local Mach number and temperature 

ratio are presented as computed in this manner. It will be noted that 
the curves exhibit a slight change for different Reynolds number which 
was due to a change in calibrated tunnel Mach number at different pressure 
levels . 


Laminar Recovery Factor 

The calculated values of the ratio H T and of the laminai- 

recovery factor are shown in figure 9 * The recovery factors were cal- 
culated by a method of Van Driest (ref. 11) and by q r = These 

values are compared with data from the recovery-temperature model. 

Included are the recovery temperature ratios and the recovery factors 
calculated from the uncorrected recovery temperatures or equilibrium 
temperatures, and also from the recovery temperatures corrected for 
conduction and radiation. 

The laminar recovery factor was also calculated from q r = Jfr, 

where the Prandtl number was calculated for the average of the total and 
the local stream temperature as was suggested by Stine and Wanlass 
(ref. 3 )» The results of this calculation agreed closely with those 
from the method of Van Driest and were not included in figure 9. 

The experimental values of T r /T t were corrected for conduction and 

radiation by first calculating the local heat flow into the model at equi- 
librium temperature and then by using this value of q s to obtain a value 

°f Plot on respective curves of d s /T t against T/T^. In 
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this way a point was obtained sufficiently near the T/Tj. axis to enable 
the fairing from the heat-transfer-model data to = °* ®i e value 

of T r ^T t obtained by this method was the recovery temperature corrected 
for radiation and conduction. 

The local heat flow q s into the model at equilibrium temperatures 
was calculated by using the following equation: 


% 


J t b ~ r in\ 5 
Ir^sin 0 J ^0 


HI 


+ ere 




(7) 


The first term on the right represents the heat flow into the model as 
found by equation (6) for a model in equilibrium, that is, = 0. The 

second term on the right represents the heat flow into the model which 
is lost by radiation. In a slightly different notation this relation is 
shown in reference 17. The temperature of the surroundings T a was 

assumed to be constant at lAO° F over the entire hemisphere. 

In order to apply equation 7, an estimate of the emissivity must be 
made. This estimate depends upon the surface condition of the material 
as well as the kind of material and its temperature. At the start of 
the testing the model had a polished surface, but after the testing the 
surface of the hemisphere had become slightly roughened. This roughness 
was greatest at the stagnation point and was slight at the sphere-cylinder 
juncture . Because of the uncertainty in the value of the emissivity at 
the elevated temperatures of these tests , it was assumed to be constant 
(0.15) over all the model. It was felt that this estimated value of 
emissivity was not in error by more than ±0.05. 

When the corrected values of temperature ratio T r ^ T Q are compared 

with the values obtained from the two theories, the experimental data 
does not decisively select either theory. Since the Van Driest method 
considers a local Mach number effect on the recovery temperature ratio, 
the values of T^ T q obtained by the method were used in the reduction 

of the heat -transfer data. 


Temperature Distribution 

Surface-temperature distributions obtained for various times during 
a typical run on the heat-transfer model are presented in figure 10. At 
the nose the average value of the temperature was about 700° R and on the 
cylinder about 600° R giving values of Tg/T*, of approximately 7 and 6, 
respectively. 
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Heat-Transfer Parameters Nu/ y^Re 

The local laminar heat-transfer parameter was calculated on the 
basis of two different characteristic lengths. For one parameter 
Nu x / /Rex, the characteristic length was the distance x along the body 
generator from the stagnation point. The air -flow properties used in 
this parameter were the local flow conditions just outside the boundary 
layer. The values of this parameter for M = 6.8 are shown in figure 11 
along with the range of experimental data at M = 1.97 from reference 3. 
These experimental values are compared with theories by Stine and Wanlass 
(ref. 3), by Chapman and Rubesin (ref. l 8 ), with a modification of a 
theory by Sibulkin (ref. 7), and with a theory by Reshotko and Cohen 
(ref. 12). For the other heat -transfer parameter Nuj/ Re-p the diam- 
eter of the model was used as the characteristic length, and the air-flow 
properties used were taken as those immediately behind the normal shock 
in the center of the bow wave. The values of this parameter for M = 6.8 
as well as the range of values for M = 1.97 (ref. 3) are shown in fig- 
ure 12. These values are compared with the theory by Stine and Wanlass 
and with a second modification of the theory by Sibulkin. 

The local heat-transfer parameter Nux/y/Rex, as calculated by the 
theory of Stine and Wanlass, is based upon the local Mach number distri- 
bution at the outer edge of the boundary layer and upon the local Prandtl 
number. In the present tests, the Prandtl number varied from about 0.68 
at the stagnation point to about 0.7^ near the end of the cylinder. A 
value of Pr = O. 7 O was used throughout the calculation of the theory 
of Stine and Wanlass. This caused an estimated maximum deviation from 
the mean heat-transfer parameter of only 2 percent. 

In the Stine and Wanlass theory, the assumption is made that the 
difference between the actual wall temperature and the insulated wall 
temperature is small. In the present tests, this temperature difference 
varied from about 100° F to about 500° F. An examination of figure 11 
discloses that the data points have a tendency to fall below the theory. 
This difference may be due in part to the rather large temperature dif- 
ference. The reasonable agreement of the data with the theory, however, 
shows that the theory of Stine and Wanlass may be used with good results 
even when the insulated wall temperature and the actual wall temperature 
are quite different. 

In comparing the theoretical results of Stine and Wanlass at M = 6.8 
and at M = 1.97, a slight discrepancy appears in that a discontinuity 
occurred at the hemisphere-cylinder juncture (x/b = 1.00) for M = 6.8 
and not at M = 1.97* This discontinuity can be traced to a difference 
in local Mach number distributions obtained from the two investigations. 

In this investigation, the pressure distribution has an abrupt change in 
slope at the juncture (from a high pressure gradient to a very low pres- 
sure gradient) which gives a similar change in the local Mach number 
gradient (fig. 7 ) • In reference 3 , no such change was noted in the 
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pressure distribution, and thus the local Mach number gradient had a 
gradual change at the juncture resulting in a continuous curve for the 
local heat-transfer parameter based on x. Theoretically, a discontinuous 
pressure gradient is predicted at the juncture. 


Although the use of the theory of Stine and Wanlass leads to the 
calculation of the local heat-transfer parameter Nu x j Re x , the param- 
eter N u p/ sjltej) may be calculated from the parameter based on x, when 


the free-stream Mach number, the local Mach number on the model, and the 
stagnation temperature are known. When the experimental and theoretical 
values of the parameter Nu in reference 3 were converted to 

Nu jRe-£)> a stagnation temperature of 520° R was assumed. 


The theory of Chapman and Rubes in (ref. 18) is a method of obtaining 
the values of the local heat-transfer parameter for a flat plate. The 
air -flow properties used were those at the outer edge of the boundary 
layer. The distance along the body generator was taken as the flat-plate 
distance from the leading edge. This theory depends directly upon the 
air temperature at the outer edge of the boundary layer and upon the 
surface temperature and is not directly a function of the Mach number. 

The theory evaluation at M = 1.97 is that shown in reference 3. 

In reference 7 , Sibulkin suggested applying his incompressible 
solution for stations near the stagnation point of a sphere to the nose 
of a blunt body moving at supersonic speeds by using as free-stream 
conditions those behind the center of the bow wave. Since the velocity 
gradient at the nose would not be expected to be the same as in incom- 
pressible flow, Sibulkin's equation should be used before the incom- 
pressible velocity gradient around the sphere is incorporated, which 
gives 


Nu l = 0.763Pr 0 , ^L ^ 0,5 


(8) 


The relation was derived with the assumption of incompressible flow, but 
the author suggests that the theory may be valid for the flow near the 
stagnation point of a blunt body in a supersonic stream. 

Equation (8) may be written as 

4 Re L 


hL 

k 



=-0.763Pr 


O.k /gL^O. 


(9) 
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When the characteristic length is x and the flow velocity and other 
physical conditions of the air are evaluated at the outer edge of the 
boundary layer, the first modification of Sibulkin's theory is obtained. 
This appears as follows: 



= 0.763PT 0,4 (10) 


The theoretical heat-transfer parameter evaluated in this manner is a 
function only of the Prandtl number. 

When the characteristic length is D and the flow velocity and other 
physical properties of the air are evaluated immediately behind the normal 
shock, equation (9) becomes 

hD 

^ U D _ ^0 

/pQ^gP 

\l da 


= 0.763Pr 


0.4 


(PD\0.5 

\ U °J 


( 11 ) 


This equation is the second modification of Sibulkin's theory used in 
this report. The theoretical heat-transfer parameter evaluated in this 

manner is a function of the Prandtl number and of — . 

u g 

In figure 13(a), the values of are plotted against x/d 

on the hemisphere. These curves show that 3 is constant near the 
stagnation point at a given Mach number. In figure 13(b), experimental 
values of 3D/u ct from several references, as well as curves computed by 
using modified Newtonian and incompressible flow theories, are plotted 
against Mach number. This curve shows that 3D/u CT is a function of 

Mach number. Therefore, Nuj^/^Re^, as found by Sibulkin's theory near 
the stagnation point, is a function of Mach number . 
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When the Newtonian theory was used to find 0 D/u ct , the Newtonian 

pressure coefficient as shown in figure 5 was used to determine the 
pressure distribution around the nose of the hemisphere. The correct 
pressure coefficient at the stagnation point to he used with each Mach 
number was found from normal -shock-flow relations (ref. 16 ) . The 
isentropic-flow relations between the pressure at the stagnation point 
and the surface pressures computed for subsequent stations were used to 
calculate the local flow conditions just outside the boundary layer. 

This is the same procedure that was used to calculate the local flow 
conditions from the measured pressures. 

In figure 11 it is seen that the data indicates a valu e at the stag- 
nation point about 10 percent below the value of Nu^J^ Re x as found by 

Sibulkin's method. In this figure it is seen that the data from tests 
at two different Mach numbers, M = 1.97 and M = 6 . 8 , agree within the 
scatter of the data. In figure 12, it is seen that the data at the free- 
stream Mach numbers of 1.97 and 6.8 indicate values wh ich, within the 
scatter, are consistent with the predictions of Nu \J^- e Y) dy Sibulkin's 

method. It is important to note, however, that the variation of 0D/u ct 

with the free-stream Mach number must be taken into account and that the 
prediction of 0 by the potential theory may not be used for a hemi- 
spherical nose subjected to supersonic flow. 


Local Stanton Number 

In figure l4 the local Stanton number is plotted against the local 
Reynolds number. The data presented includes points from 0=5° to 
0 = 4 5 ° on -the hemisphere surface. 

By assuming that the velocity and temperature boundary-layer profiles 
were similar and of parabolic shape, Sibulkin (presented in ref. 2) derived 
an approximate solution for Nu D over the entire hemisphere for incom- 
pressible flow which, when transformed for the local conditions outside 
the boundary layer and distance from the stagnation point, becomes 


Nu = k(9)x^j^ ^Pr 


0.4 


( 12 ) 


where k( 0 ) is a proportionality function which depends on the angular 
position of the point under consideration measured from the axis of flow 
symmetry . 

If 0 is assumed constant (a valid assumption for 0 up to 45° in 
these tests, see fig. 13 ), the equation reduces to 
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(13) 


From reference 6 it may be seen that k(9) may be approximated near the 
stagnation point by the value of O. 763 , and the local Stanton number 
near the stagnation point becomes 


For this equation to give the correct Stanton number around the hemisphere 
surface, du/dx or 3 must be constant and k(0) should not depart 
radically from the value of 0.763* 

In figure 14 the heat-transfer data has been compared with Sibulkin's 
stagnation-point solution shown in equation (l4 ) . As seen from fig- 
ure 13(a), p can be considered constant up to values of 9 greater 


k(9) may be obtained from the decrease of the heat-transfer parameter 
from the constant stagnation -point value as x is increased. Indeed, 
even the difference in level which is observed between the stagnation- 
point curve and the data in figure 14 may also be detected in figure 11. 
Sibulkin's stagnation -point value is approximately 9 percent above the 
data at the nose and about 15 percent above the data at the 45° station. 
Above 45°, k(9) diverges rapidly and 3 begins to diverge. Inasmuch as 

the data above 45° do not correlate, they are not presented in figure 14. 


The average heat-transfer coefficient determined experimentally in 
this investigation and by Stalder and Nielsen in reference 4 are presented 
in figure 15. The data from Stalder and Nielsen are those for 
M = 0.12 to 0.17, M = 1.75, and M = 2.67. Also included in figure 15 
are experimental points found by an integration of experimental data by 
Stine and Wanlass (ref. 3) at M = 1.97* The data has been compared with 
values found by using the theory of Stine and Wanlass at M = 1.97 and 
M = 6.8 and with the incompressible theory of Sibulkin used by Korobkin 
(ref. 2). Korobkin adjusted the latter solution to pass through the 
stagnation point for incompressible flow (ref. 7)* 




(l4) 


than 



In figure 11 an indication of the variation of 


Average Heat -Transfer Coefficient 
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The data from Stalder and Nielsen (ref. A) for Mach numbers below 
2.67 appear to correlate with the adjusted incompressible theory. The 
data of reference A were evaluated, however, in a different manner from 
the data of this paper. In the present paper the local values of the 
heat-transfer coefficient were calculated from the local heat convection 
and the local recovery temperature. Then these local values of h were 
integrated over the area of the hemisphere to find the average value. 

In reference 4 the average heat-transfer coefficients were evaluated by ^ 
measuring the total heat convection to the hemisphere and an "effective" 
average recovery temperature found by an extrapolation to zero net heat 
transfer to the surface. The different results as found by the two 
methods are seen in figure 15- 

The integrated data from Stine and Wanlass and the data from this 
investigation show a good correlation at Mach numbers of 1.97 and 6.80 

and Reynolds numbers from 0.002 x 10^ to 2.28 X 10^. 


Schlieren Observations 

The normal shock and body boundary layer are shown in the schlieren 
photographs (fig. l 6 ) for a single-pass and double-pass schlieren system. 
The relatively nonsensitive single-pass system clearly shows the normal 
shock while the sensitive double-pass system was used to show that the 
boundary layer is laminar on the surface of the cylinder. The two flow 
photographs appear to show additional disturbances about three-fifths of 
the distance along the cylinder from the hemisphere -cylinder juncture. 
These disturbances result from the intersection of the bow shock wave 
and the boundary layer at the windows . The heating effect was great 
enough so that a similar picture of this region was obtained immediately 
after a run, as shown in the post-run zero picture for the double-pass 
system . 


CONCLUSIONS 


A program to investigate the aerodynamic heat transfer of a 
noni sothermal hemisphere -cylinder has been conducted in the Langley 
11-inch hypersonic tunnel at a Mach number of 6.8 and a Reynolds number 

from approximately 1.09 x 105 to 1.03 X 10 6 based on diameter and free- 
stream conditions. In this investigation the boundary layer was laminar 
having average values of T s / T^ at the nose of about 7 about 6 on 
the cylinder (where T s is the local free-stream temperature just outside 
the boundary layer and Too is the free-stream temperature ahead of the 
normal shock) . The data have been correlated with theoretical analyses 
and results of other experiments at lower Mach numbers. The following 
results were obtained: 
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1. The experimental heat-transfer coefficients from this investi- 
gation on a nonisothermal (constant wall thickness) hemisphere -cylinder 
with heat flowing to the surface were slightly less over the whole body 
than those predicted by the theory of Stine and Wanlass (NACA Technical 
Note 33^) for an isothermal surface. 

2. A modification of Sibulkin's stagnation point solution gave the 
trend of the local Stanton number with local Reynolds number for angles 
up to 45° from the stagnation point. The calculated values., however, 
were approximately 12 percent higher than the experimental values . 

3. Pitot profiles taken at a Mach number of 6.8 on a hemisphere- 
cylinder have verified that the local Mach number or velocity outside 
the boundary layer required in the use of the theories may be computed 
from the surface pressures by using isentropic flow relations and condi- 
tions behind a normal shock. The experimental pressure distribution at 
a Mach number of 6.8 is closely predicted by the modified Newtonian 
theory. The velocity gradients calculated by using modified Newtonian 
theory at the stagnation point vary with Mach number and are in good 
agreement with those obtained from measured pressures for Mach numbers 
from 1.2 to 6.8. 

4. At the stagnation point, a second modification of the theory 
by Sibulkin using the diameter and conditions behind the normal shock 
was in good agreement with the experiment when the velocity gradient at 
the stagnation point appropriate to the free-stream Mach number was used. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 21, 1956 
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Figure 1.- Typical construction of models 
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(a) Heat-transfer and pressure models. 
Figure 2.- Hemisphere -cylinder models. 
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(b) Recovery-temperature model with 2-inch diameter and 0.005-inch-thick wall. 

Figure 2.- Concluded. 
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Figure 3 *- Heat transfer to the hemisphere of the 3 •025-inch-diameter 
heat-transfer model both with and without the conduction correction. 
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Figure A.- Heat -transfer rate plotted against surface temperature 
various angles of 0 on the 3- 025-inch-diameter heat-transfer 


for 

model. 
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Figure 5.- Longitudinal pressure-coefficient distribution. 
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Figure 6.- Pitot surveys of flow near the surface of the cylinder 
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Figure 8.- Variation of the local temperature ratio (outside the boundary 

layer) along the surface. 
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Figure 9*- Variation of the surface temperature ratio and recovery factor 

along the surface. 
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Figure 10,- Typical surface temperature distribution on 3«025-inch-diameter 
heat-transfer model at various times during a run. 
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Figure 11.- Variation of the local heat -transfer parameter along the surface. 
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Figure 12.- Variation of the heat-transfer parameter (based on diameter and conditions behind 

the normal shock) along the surface . 
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(a) Nondimens ional local flow velocity along the surface. 
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(b) Nondimens ional velocity gradient at the stagnation point of the hemi- 
sphere plotted against free-stream Mach number. 

Figure 13.- Nondimens ional flow parameter on the hemisphere outside the 

boundary layer. 
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Figure 14.- Variation of the local Stanton number with the local Reynolds number for values of 9 

from 0° to 45°. 
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Figure 15 •- Average heat transfer on the hemisphere plotted against Reynolds number for conditions 

behind a normal shock. 
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Figure l6.~ Schlieren photographs 


of shock and boundary layer on hemisphere-cylinder at 
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